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The spectrum of two-gluon glueballs below 3 GeV is inves- 
tigated in a potential model with dynamical gluon mass using 
variational method. The short distance potential is approxi- 
mated by one-gluon exchange, while the long distance part is 
taken as a breakable string. The mass and size of the radial 
as well as orbital excitations up to principle quantum number 
n = 3 are evaluated. The predicted mass ratios are compared 
with experimental and lattice results. 

PACS numbers: 12.39. MK 



Quantum Chromodynamics (QCD) is widely accepted 
as the theory of strong interactions. It is generally be- 
lieved that the gluon self-coupling in QCD implies the 
existence of bound states of confined gluons known as 
glueballs. The experimental discovery of these glueballs 
would be very important and would give further support 
to the theory of QCD. However, numerous technical dif- 
ficulties have so far hampered our unequivocal identifica- 
tion of glueballs by experiment, largely because glueball 
states can mix strongly with nearby qq resonances. Nev- 
ertheless, the estimation of mass and size of pure gluon 
glueball states should still be pursued. This could guide 
experimental searches, as well as provide calibration for 
models of glueballs. 

Over the past 20 years there has been an on-going effort 
to obtain a nonperturbative form for the gluon propaga- 
tor. Perhaps one of the most interesting result is that 
the gluon may have a dynamically generated mass |Q. 
The existence of a mass scale, or the absence of a pole at 
k 2 = 0, is natural if one assumes that gluons do not prop- 
agate to infinity, i.e., these propagators describe confined 
gluons. The concept of massive gluon has been widely 
used in independent field theoretic studies, and examples 
about the consequences of massive gluons can be found 
in the literature . 

In this paper, we focus on the calculation of two-gluon 
glueball systems and extend our previous work S on 
the estimation of the mass and size of low lying glue- 
ball states, using the variational method in the potential 
model of Cornwall and Soni P,^0[. The main feature of 
the present work is the consideration of radial as well 
as orbital excitations, up to principle quantum number 
n = 3. 

To exhibit both asymptotic freedom and the non- 
Abelian nature of QCD, gluon dynamics can be described 
as massive spin-one fields interacting through one-gluon 



exchange and a breakable string. At short distance the 
effective coupling constant of the gluon-gluon interaction 
becomes small and the interaction can be treated pertur- 
batively. The short distance potential is approximated 
by one-gluon exchange and can be extracted from the 
tree-level Feynman amplitude of Fig. 1, 
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where q is the momentum transfer of the system. At long 
distance, the non-Abelian nature of QCD implies gluon 
confinement via nonperturbative effects. These nonper- 
turbative effects are implemented by introducing a string 
potential V^tr which is assumed to be spin independent, 

VW = 2m(l - e-? mr ), (2) 

where (3 is related to the adjoint string tension Ka via 
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In the potential V s tr, the color screening of gluons is 
brought about by a breakable string; that is, the adjoint 
string breaks when sufficient energy has been stored in 
it to materialize a gluon pair. This form of the string 
potential simulates the inter-gluonic potential as seen in 
lattice calculations [ pr| . 

Thus the gluon-gluon potential relevant for two-gluon 
glueballs is M 
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where A is defined as 
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and is related to the strong coupling strength of the pro- 
cess. The terms containing A are from the short distance 
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FIG. 1. Diagrams for gg — * gg scattering. 
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potential. Note that the color wavefunctions have been 
contracted out with the structure constants in gluon ver- 
tices [|| , so that the strength of the gluon-gluon coupling 
is three times as large as that of the gluon-quark cou- 
pling. In Eq. (0), S = Si + S 2 is the total spin of the 
2-gluon glueball. Note also that Si acts on the polariza- 
tion vectors ei and e3 while S2 acts on the polarization 
vectors e 2 and e4. Each spin operator S„ = (S 1 , S 2 , S 3 ) 



is defined as (S )i 



-ie lJ , which satisfy 



[S^S^ = ie ijk S k 
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i.e. S 1 , S 2 and S* 3 are SU(2) group generators as desired. 

For this gluon-gluon potential we are left with three 
parameters: effective gluon mass m, string breaking pa- 
rameter (3 and adjoint strong coupling constant A. In this 
model the constituent gluon mass is evaluated to be 600- 
700 MeV, which is roughly twice the constituent quark 
mass p| . In the intermediate or at least at long distances, 
the fundamental string tension Kp is one of the most fun- 
damental physical quantities in quark confinement, and 
is related to the Regge slope a' by 12 



K F = 
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= 0.18 GeV 2 
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with the experimental value for a' ~0.9 GeV -2 . The 
adjoint string tension Ka for gluon confinement can be 
related by the strong evidence of Casimir scaling hypoth- 
esis lOl on the lattice via 
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Thus, the string breaking parameter (3 is about 0.4-0.6. 
For the gluon propagator with a dynamical mass the ad- 
joint strong coupling constant A at one-loop level turns 
out to be fjj 
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where £ ks 4. The higher order correction at two loop 
level does affect the value, as shown in the case of the 
fundamental strong coupling constant a s (Q 2 ) How- 
ever, it has not been solved so far in the literature. Nev- 
ertheless, A is expected to be in the range of 1.5 ± 0.5. 

The invariance of charge conjugation and parity for 
strong interactions implies that a two-gluon glueball must 
have quantum numbers J PC = J^ -1 ) ( + '. On the other 
hand, since the gluon is a spin 1, color octet boson, the 
color singlet wavefunction of a 2-gluon glueball is sym- 
metric. Hence, a symmetric spin wavefunction with spin 

or 2 must be accompanied by a symmetric spatial wave 
function with an even value of orbital angular momen- 
tum, and an antisymmetric spin wavefunction with spin 

1 must be accompanied by an antisymmetric spatial wave 
function with an odd value of orbital angular momentum. 
In this way one can count the number of states. For ex- 
ample, we have two S-wave, three P-wave and six D-wave 



states with definite J PC corresponding to principle quan- 
tum number n — 3. 

The Hamiltonian of two-gluon glueball system is 
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H = 2m- — V 2 + V 2g . 
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Since the color wavefunction has been contracted out, the 
remaining wavefunction is of the form 

^nljm(r) = \nljm) oc t/j n i(r)\jm) 
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where 
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The trial radial wavefuctions ipni( r ) are constructed by 
orthogonality as follows, 
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where a is the variational parameter. The spin wave- 
function Xsm 3 is constructed by the direct product of 
two gluon polarization vectors e^ Al ' and e^ 2 '' and, in 
turn, the total angular momentum eigenstate \jm) can 
be constructed by the direct product of orbital eigen- 
state, the spherical harmonics Yi mi (9, (f>), and the spin 
wavefunction Xsm 3 - The coefficients (lAilA2|sm s ) and 
{lmism s \jm) in Eqs. ( pi] ) and ( |l2] ) are just the Clebsch- 
Gordan coefficients. 

When considering the case with L = JH], we drop 
the spin-orbit and tensor terms in Eqs. (ffl^ However, 
they do contribute for the case with L > 0. Since only 
the tensor term in Eq. (^) is related to the total angular 
momentum eigenstates \jm), we need to calculate their 
expectation values. Defining the tensor operator as, 



T = S 2 -3(S-f) 2 



(14) 



integrating out the spherical harmonics and doing spin 
algebraic calculation on spin wavefunctions, we obtain 



(0 0|T|0 0) = 0, for 1 = 0, s = 0, 
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and 



(2 to|T|2 m) = 0, for 1 = 0, s = 2. (16) 
For I — 1 and s — 1 we have 
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(2m|T|2m) = -, 
5 

(lm|T|lm) = -1, 
(0 0|T|0 0) = 2. 

For I — 2 and s — 0, we have 

(2 m|T|2 m) = 0, 

and for / = 2 and s — 2, we have 
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Hence the contribution to the glueball mass M can be 
written as 



M/m = 2 + E K + E Y 
where 

1 



E 5 +E LS + E T + E stI , (20) 
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and \nl) = ip n i(r). 

We then use the variational method with trial radial 
wavefunction ip n i (r) to evaluate the glueball mass M and 
root-mean-square radius r rms = yj (r 2 ) for each glueball 
state. The lightest scalar and tensor glueballs have been 
investigated in Ref. For the lightest scalar glueball 
there is an attractive <5-function term in the potential and 
hence the Hamiltonian is unbounded from below. It has 
been conjectured that this maximum attraction channel 
in ++ state could be related to the gluon condensation 
that triggers confinement. We proposed a physical solu- 
tion by smearing the gluon fields; that is, we replace the 
5 function by the smearing function 



which approaches (5 3 (r) for k — > oo. In contrast to the 
lightest scalar glueball, the lightest tensor glueball is sta- 
ble since the (5-function term is repulsive. As the vari- 
ation is slight in the mass and size estimation (^] of a 
glueball in the A-/3 parameter space, we take the central 
values of A =1.5, and (3 =0.5 fllq] . 

The lattice result of M^Sq), M{1 5 S 2 ) = 1730, 2400 
MeV is taken as input for the lightest ++ and 2 ++ glue- 
balls, respectively From the 2 ++ input, we find 
the constituent gluon mass ~ 670 MeV, about twice the 
constituent quark mass, and the lightest tensor glueball 
is found to have the typical hadron size of ~ 0.8 fm. 

The value for k in Eq. (E3) is fixed by the mass ratio 
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1.39, 



(23) 



from the converging lattice and experimental results. 
The size of the lightest scalar glueball is found to be a 
mere ~ 0.1 fm. We note that, although we always have 
an attractive (5-function term in the potential for scalars, 
only for the lightest scalar glueball is the smearing of 
gluon fields needed. For all other scalars, the (5-function 
potential gives very small mass corrections. One can con- 
sider the sum of kinetic energy and (5-function terms as 
the effective kinetic energy. We find that the sum con- 
tributes less than 12% to all glueball masses, except for 
the lightest scalar glueball, which is 43%. This may be 
marginal for a nonrclativistic treatment, and stretching 
the applicability of our relativistic expansion. 

We list the mass spectrum and size of two-gluon glue- 
balls up to n = 3 in Table 1. 

Table 1. Masses and sizes of two-gluon glueballs up to 
n = 3, with the lightest ++ and 2 ++ masses taken as 
input. 
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0++ 


{1730} 


0.1 
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2++ 


{2400} 


0.8 
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2++ 


2700 


1.6 
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2690 
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2++ 


2693 
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3++ 


2694 
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4 ++ 


2695 


1.7 
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From Table 1 we see that the glueball mass M de- 
creases with increasing orbital angular momentum L for 
fixed n. For fixed n and L = 0, 1 (2), M increases (de- 
creases) with total spin angular momentum S. In turn, 
both M and size increase with total angular momentum 
J for fixed n, L, and S. On the other hand, M increases 
with n at fixed L and S. All these increments are slight, 
except for the n = 1 case where the attractive S potential 
is present. 

Our calculations show that glueball masses are almost 
independent of A, and increases with (3 only slightly |l5"| . 
Although our pure gluon glueballs are different from real 
glueball states which can mix strongly with nearby qq 
resonances, the better way to compare with experimental 
data is to take mass ratios, either to eliminate the A and (3 
dependence, or to reduce the effect of mixing. We list in 
Table 2 the glueball masses in increasing order, together 
with their corresponding mass ratios with respect to the 
lightest tensor glueball. Comparison with Table 1 shows 
that, for a given quantum number, the ordering in mass 
is also the ordering in size. 

Table 2. Comparison of predicted glueball masses (mass 
ratios, normalized to lightest 2 ++ ) with lattice |L6| and 
sample experimental data. The superscripts a, b, c and 
d indicate data coming from ITS, pll, Jl9] and p3], re- 
spectively. The experimental numbers are not meant to 
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are taken as 


input. 
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MeV units. 




J PC 


Constituent 


Lattice 


Experiment 




{1730} 


(0.72) 


1730 (0.72) 


1500" (0.76) 


0++ 


2685 


(1.12) 


2670 (1.11) 


2105 ft (1.06) 




2710 


(1.13) 




2320 c (1.17) 




2790 


(1.16) 








{2400} 


(1.00) 


2400 (1.00) 


1980" (1.00) 




2693 


(1.12) 


3290 (1.37) 


2020 d (1.02) 


2++ 


2700 


(1.13) 




2240 d (1.13) 




2730 


(1.14) 




2370 d (1.20) 




2810 


(1.17) 
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(1.15) 
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(1.09) 








2770 


(1.15) 






2- + 


2615 


(1.09) 


3100 (1.29) 


2040" (1.03) 




2775 


(1.16) 


3890 (1.62) 


2300 d (1.16) 


1++ 


2690 


(1.12) 




2340 d (1.18) 


3++ 


2694 


(1.12) 


3690 (1.54) 


2000 d (1.01) 










2280 d (1.15) 


4 ++ 


2695 


(1.12) 


3650 a (1.52) 


2044 d (1.03) 










2320 d (1.17) 
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(1.58) 


3640 (1.52) 






) 3680 


(1.53) 


3850 (1.60) 




3fl(3— 


) 3690 


(1.54) 


4130 (1.72) 





Although our predictions rely on the inputs of the 
lightest scalar and tensor glueball masses, 1730 and 2400 
MeV respectively from lattice calculation, we find that 
our predicted mass ratios do find experimental correspon- 
dence, sometimes even better than comparing with lat- 
tice calculations. This may be due to the rather straight- 
forward physical picture of a potential model. However, 
the experimental situation is far from settled, and com- 
parison with lattice is quite necessary. 

As stated, the lightest tensor glueball is much more 
stable within the model than the lightest scalar glueball, 
hence all mass ratios are normalized to this state. We find 
that, except for the peculiar lightest ++ state, all ratios 
of two-gluon glueball masses are of order 1, and mass 
ratios of three-gluon to two-gluon glueballs are of order 
3/2, respectively. This is consistent with a constituent 
picture. We then need to understand the difference from 
lattice spectrum, which suggests a considerably higher 
excitation energy. 

Going back to Table 1, we note that the con- 
stituent picture gives certain multiplet structures gov- 
erned mostly by n and L, with relatively small splittings 
in S (except for n — 1). We note that there are spin-orbit 
and tensor forces at work. Thus, the "first excitation" 
from the n = 1 states seem to be the (n, L, S) = (2, 1, 1) 
(rather than (2, 0, 0)) states of 0"+, 1"+ and 2 _+ . The 
masses are all around 2600 MeV, or ~ 1.08 times the 
lightest 2 ++ . They are also of similar size as the lightest 
2 ++ , the h actually smaller, hence they are likely more 
relevant in production processes. 

The h may be most interesting since it is the light- 
est pseudoscalar glueball. Its mass of 2570 MeV is in re- 
markable agreement with the lattice result of 2590 MeV. 
The heaviness suggests [^lj that there may be little glue- 
ball mixture in rj and 77', and OZI rule violation in these 
mesons arise from vacuum effects. The 1 h state with ex- 
otic quantum numbers cannot be a qq meson and is also 
very interesting. Again, its heaviness compared to hy- 
brid candidates such as /5(1600) |^2| suggest that search 
for lightest hybrids are not complicated by presence of 
glueballs. For the 2"+, the mass of 2615 MeV in con- 
stituent model is in contrast with the lattice result of 
3100 MeV, which has a 500 MeV excitation energy with 
respect to the lowest lying h . This large gap is absent 
in our constituent model, which can largely be traced to 
the tensor force. 

The (n, L) = (2, 0) and (3, 0) states of 0++ and 2++, 
are straightforward radial excitations of the n = 1 states. 
Since the lightest ++ is especially light, the excitation 
energy is more than 900 MeV. But for the 2 ++ , it seems 
that the radial excitation energy is only 300 MeV or so, 
while going on to n — 3 states, the excitation energy 
is less than 100 MeV. We notice that the size has also 
reached beyond 2 fm, hence the drop in excitation energy 
reflects the approach to string breaking beyond a couple 
of fermis. Such behavior, however, is not seen on the 
lattice, where the second 2 ++ state is also almost 900 
MeV higher than the lightest one, which again appears 
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more "stringy" ("Regge"). 

The (n, L, S) = (3, 2, 0) and (3, 2, 2) states in 
the constituent model are close to degenerate, which 
arise again due to a balance between spin-orbit and ten- 
sor forces. Interestingly, the ++ constructed out of 
L = S = 2, possible only starting with n = 3, is slightly 
lighter than the radial excitation of the lightest ++ . It 
is thus actually the second lightest such state, with size 
of 1.1 fm, comparable to the lightest 2 ++ but only half 
that of the radially excited ++ . We note that its mass of 
2685 MeV is in excellent agreement with the second ++ 
state from lattice, although our radial excitation states 
are also not in disagreement. The (n, L, S) = (3, 2, 0) 
2 ++ state is slightly lighter and smaller in size as well 
than the radially excited 2 ++ . The interpretation of ex- 
citation energy for ++ versus 2 ++ therefore could be 
rather different between our model and lattice calcula- 
tions. 

The above two special L — 2 states are accompanied by 
the host of 1 ++ , 2 ++ , 3 ++ and 4 ++ glueballs which are 
comparable in mass and gradually growing in size. They 
could be more interesting than the "usual" radial excita- 
tions of the lightest 0++ and 2++, or 0"+, 1"+ and 2~+, 
which are two and a half fermi in size and rather large. 
In contrast, the lattice 3 ++ and 4 ++ states are close to 
3700 MeV, again appearing as stringy excitations. 

The experimental results are rather uncertain at 
present, but there does seem to be many states in the 
2000 MeV region as compared to lattice suggestions. 
In part for reasons of comparison, we take here 1980 
MeV for a 2 ++ candidate mass |L8| . Glueball candidates 
abound but they are very hard to pin down. In contrast 
to lattice results, especially for the glueballs with size less 
than 2 fm, the richness of our spectrum should give some 
hope for experimental search, although one clearly has a 
long way to go. For the lattice calculations, it is not clear 
to us whether the short distance spin-orbit and tensor 
interactions, of great importance to our model (though 
mysteriously balancing each other), are replicated. 

In our model, the size of the lightest scalar glueball is of 
order 0.1 fm. The extreme smallness may be an artifact of 
the treatment of the attractive 6 potential, but we do ex- 
pect a physically smaller lightest scalar from the heuristic 
point of view, because of the extra attraction. A more di- 
rect calculation of the ++ glueball mass and size on the 
lattice would require relatively fine lattice spacings p3[ . 
It would be interesting to see if our result of small (v + 
size could be born out on the lattice. The lightest pseu- 
doscalar glueball mass is 2570 MeV. However, there are 
the other two h states with mass 2765 and 3780 MeV 
composed of two- and three-gluons, respectively. Mix- 
ing among them should make the lightest ^ two-gluon 
glueball mass lighter. On the other hand, despite the lat- 
tice QCD prediction that the lowest lying 1 h glueball 
has mass heavier than J/-0 p4j , one should not exclude 
the 1"+ glueball in search below 3 GeV The 1~ 

and 3 only exist in three-gluon glueball states, hence 
their masses should be heavier than 3 GeV. Scaling from 



the lattice result of M 2+ + =2000-2400 MeV |16|j2J,|2j the 
mass range of these glueballs is 3.1-3.7 GeV, right in the 
ballpark of J/tjj and masses. The proximity of the 1 _ ~ 
glueball to J/tp or perhaps ip 1 may be what is needed from 
comparison of J/ip and tp' two body hadronic decays [^8) . 

In conclusion, two-gluon glueballs have been studied in 
a potential model with constituent gluons. The potential 
is approximated by one-gluon-exchange plus a breakable 
string. The mass and size of radial and orbital excita- 
tions, up to principle quantum number n — 3, are eval- 
uated by variational method. With the lowest lying ++ 
and 2 ++ masses from lattice taken as input, all masses 
of two-gluon glueballs are found to be below 3 GeV with 
size less than 3 fm. The predicted masses for the second 
++ as well as the lowest h are in excellent agreement 
with lattice. Further excitation patterns, however, differ 
considerably. While lattice calculations find fewer states 
populating 3000 MeV range, hence an excitation energy 
of 500-900 MeV, our constituent model gives two more 
++ states, four more 2 ++ states, and an additional h 
state, which are 12%-17% heavier than the lightest ten- 
sor. This pattern may well be more consistent with ex- 
periment. There are also a number of 1 =F+ , 2 h , 3 ++ and 
4 ++ glueballs. Besides a rather small size for the lightest 
scalar, there is one ++ , 2 ++ , 1 h and 2 h each 
that are 1 fm in size, each one the lightest member for 
the given quantum number, with 1 ++ only slightly larger. 
The number of glueballs clustering at 13% or more heav- 
ier than the lightest 2 ++ are typically 2 fm in size. Mass 
and size seem positively correlated. While the model cer- 
tainly has its limitations, its heuristic nature may provide 
some help to the long quest for uncovering glueball states 
in Nature. 
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